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Abstract 

Dislocations and disclinations in a lattice of Josephson junctions will affect 
the dynamics of vortex excitations within the array. These defects effectively 
distort the space in which the excitations move and interact. We calculate 
the interaction energy between such defects and excitations, determine vortex 
trajectories in twisted lattices, and discuss the consequences for experiment. 
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I. INTRODUCTION 



A host of interesting experiments have already been performed on regular arrays of 
Josephson junctions. Tuning the ratio of the pair charging energy to the Josephson coupling 
energy in such experiments allows one to pass from a regime where the dynamics is dom- 
inated by vortex excitations, to one where it is dominated by charge soliton excitations.0 
In the vortex dominated regime one can examine such issues as the ballistic propagation of 
vorticesji! the Kosterlitz-Thouless (KT) vortex unbinding transition,^! and how the inter- 
ference of vortices is influenced by electric charges in the Aharanov-Casher effect .i In the 
charge dominated regime the parallel issues of the propagation of solitonsjl! the possibility 
of the KT charge unbinding transition! and how the interference of charges is influenced by 
magnetic fields in the Aharonov-Bohm effectffl have been investigated. The charges and 
vortices in a superconductor are approximately dual, and the physics of the two reflect this 
symmetry.@'0 

In this paper we study the motion of vortices in arrays that possess topological defects. 
Studies of Josephson junction arrays usually focus on either square or trangular lattices, 
where the structure is uniform. It is an unspoken assumption that the lattice of junctions 
be free from spatial defects such as dislocations (e.g. a missing line of junctions as in 
fig. QT].a) or disclinations (e.g. a missing wedge of junctions, as in fig.([I].b)). We focus on the 
regime where Ej, is larger or on the order of Ec, so that vortices can be viewed as stable 
massive particles with a mass determined by their charging energy. Under these conditions 
vortices are free to move ballistically in the array. We find that the defects act as distortions 
in the 2D world in which the excitations move, with two effects. First, "straight lines" 
(geodesies) in this twisted coordinate system differ from those in Cartesian coordinates, 
and second the fields surrounding the excitations are perturbed by the twisted geometry, 
alterring the interaction between vortices and leading to an interaction between vortices and 
the geometrical defects. 

The fundamental origin of the interaction between vortices and topological defects is not 
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hard to see. The self-energy of the vortex depends upon spatial gradients in the supercon- 
ducting phase. Dislocations and disclinations involve either the removal or insertion of extra 
junctions in the pristine lattice. From the point of view of the vortices this means there 
are either missing or added regions of space in the two dimensional plane. This affects the 
gradients and thus the self energy. 

In section |1| below we develop a continuum description of the dynamics on a distorted 
lattice. We then analyze the effect of topological irregularities on the dynamics of a single 
vortex in an array in which Ej > Ec, i.e., an array in which a vortex can be viewed as a 



stable massive particle. In section III we calculate the interaction of vortices or charges on 



the twisted lattice. We use these results in section [TV] to calculate the motion of vortices 
in a lattice with defects. The effects of the defects should be observable in experiments in 
which the ballistic motion of vortices is probed. We conclude in section [V|. 

In the calculations below we use Einstein summation notation in which we sum over 
repeated indices. Roman indices range over the values 1 and 2, greek indices for the 2+1 
dimensional space-time range over 0, 1 and 2. Bold face variables are 2-vectors; hatted 
vectors are normalized. 



II. DEVELOPMENT OF THE CONTINUUM DESCRIPTION 

We consider an array of Josephson junctions in an externally applied magnetic field 
perpendicular to the array. We assume that the array has only smooth distortions in its 
structure, and a small number of defects such as the dislocation or disclination in figfj]. 
Arrays of Josephson junctions are customarily described in terms of the set of variables ip a , 
and n a , the phase of the super-conducting order parameter and the number of excess Cooper 
pairs on super-conducting island a. Although the islands themselves are extended objects, 
we consider the phase and charge as being defined at a discrete set of points. (That is, we 
require the island width to be smaller than or on the order of the correlation length, so that 
the phase across the island is a constant). The thermodynamics of the array is derivable 
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from a Hamiltonian made out of two contributions, the Coulomb charging energy and the 
Josephson coupling energy. The partition function can be written as a path integral over 
the charge and phase degrees of freedom: 

Z = Y^\ / D(p a exp^--S[n a , <£> Q ]J (1) 



with the action given by, 

(2) 



S[n a , ifa] = / dr 



o 



-ihn a Lp a - 2e 2 n a n b C~ b l + E j (1 - cos(y2 a - ip b - ab )) 



o,6 

where repeated indices are summed and the path integral in r is computed in the standard 
step-wise fashion. The quantity C~\ is an element of the inverse capacitance matrix element 
between islands a and b, Ej is the Josephson coupling energy between two nearest neighbor 
junctions, and (3 = 1/ksT where T is the temperature and Icb is Boltzman's constant. The 
primed sum indicates that only nearest-neighbor phase differences are included. The phase 
difference 9 ab is caused by the externally applied magnetic field and is determined from the 
A, the magnetic vector potential, via 

6 ab = / 6 A(x)-rfx. (3) 

J a 

where we have absorbed a factor of 2e/c into A(x) 

We next make the standard Villain approximation to decompatify the phase.0 We in- 
troduce a set integer variables, defined on each island, and v ab , associated with each 
nearest-neighbor pair of islands and approximate eq.(Q) as: 

Z w c J Vn a J V(f a J2 ex P j ~ \ J Q -in a (?i(p a + 2tt^ 0) ) - 2e 2 n a n b C~ b 



(0) 



+ id i^a ~ ~ Oab + 2TrV ab ) 2 
a.b 



(4) 



where Co is a constant prefactor whose value depends upon Ej and the capacitance. The 
variable n a is now continuous; by summing over we recover only its integer contribu- 
tions. Following the Villain transformation, the phase (p ab are all regular {i.e. non-compact) 
variables. The vorticity is all contained in the variables v ab . 



The system described by eq(||), or eq(f|), is evaluated on a lattice of points, a lattice that 
may possess topological defects such as dislocations or disclinations. As in the case of a 
regular lattice, the long wavelength low frequency properties are obtained from a continuum 
description of the lattice.SS In the case of a distorted lattice, however, the continuum 
description should reflect the lattice distortions. Our goal is the construction of a continuum 
Lagrangian whose long wave length low frequency properties are the same as those of the 
discrete model of eq(£D . 

Although the phase ip a is only defined at a discrete set of points, we may define a 
continous function <£>(x) which smoothly interpolates between the phases defined at each 
lattice point. We replace the phase difference between two islands by a Taylor expansion: 

« V^(Xo) • u (i) (x a ) (5) 

where U(j)(x a ) is a vector connecting two adjacent islands. There will be two such vectors (i = 
1, 2) associated with each island, as shown in fig.£| The direction and length of the vectors 
will change as we move through the distorted lattice. In a similar fashion we associate a pair 
of the integers, v a b and v a y with these directions so that V(i)(x a ) = t> afe U(i)(x a )/w(i)(x a ) and 
V(2)(Xo) = v ab > U( 2 )(x )/ii(2)(x a ) where w w (x a ) = |u (i) (x )|, and U(i)(x a ) = u (i) (x a )/w w (x a ) 
Similarly, we approximate the contribution from the magnetic field by: 

9 ab « A(x ) • U(i)(x ) Q M « A(x a ) ■ u (2 )(x a ) (6) 

so that our partition function becomes 

Z « c / Pn(x ) f Vipi^a) ex P I ~\ I [-in(xa)(^(xa) + 2ttv^)- 

v (l) (x„),»(»)(x n ) I nJ ° 



E; 



2e 2 n(x a ) n(x 6 ) C'^ + E E [(<^( x a) - ^(x a ) + 27rof i) (x)) vfa(x) 



a i=l,2 



(7) 



We would like to turn this discrete sum over the lattice points x a and nearest neighbor 
vectors U(j)(x) into an integral over the Cartesian coordinates x\ and X2 in the plane. 



All geometrical information about the lattice is contained within the vectors U(j)(x), 
since they indicate the direction and distance to the nearest neighbors. Locally each of the 
U(j)(x) can be written as 

u (i) = (a\ ei + a\ e 2 ) (8) 

where we have suppressed the spatial dependence and the are the Cartesian unit vectors. 
We note that the a 1 • are the 2D analogs of basis triads in differential geometry; they relate 
differentials in the local, twisted reference frame (the U(j)) to the Cartesian frame.0 We may 
rewrite gradient-squared terms as: 

Yl ( V( ^ " u «) 2 = d j^ d k^P a- a* (9) 



i=l,2 

We introduce the matrix gi = a/Oj ; its inverse, denoted by gjk, is the metric tensor. A 
small displacement in the Cartesian coordinates, dx = dx 1 e± + dx 2 £2, has a distance when 
measured in the number of junction interfaces crossed (or dimensionless junction "hops") of 

ds 2 = g jk dx j dx k (10) 

whence the name "metric" . 

Furthermore, the area a of the plaquet bounded by um and U( 2 ) is just 

A = |U(i) X U( 2 )| (11) 

^et g~j u ( 2) (12) 
= g 1/2 M( 2 ) (13) 

Note that the metric tensor and its determinant may vary in space. 

We must also convert the capacitance terms to a continuum form. The original capaci- 
tance matrix satisfies Q a = C a bVb where Vi, is voltage on island 6; if we assume only nearest 
neighbor coupling, then this can be written as:0 

(4Cx + C )V a -C x Y, yb = 2en a (14) 

neighbors 



where the sites b are neighbors of a. The variables Co and C± are the diagonal and off- 
diagonal elements of the capacitance matrix. If we neglect the diagonal self-capacitance, 
(Co = 0), it states that the discrete divergence of V equals the charge on a given island. 
Through an analysis similar to what was done for the phase, we can define a continuous 
field V(x) and replace the above discrete equation by 

(d g-^d^g^d, + C (x)) V(x) = 2e n(x) (15) 

where n(x) consists of a sum of (^-functions, and Co(x) is a capacitance density given by 
Co/a. If we neglect this self-capacitance term, then the left hand side is the Laplace- 
Beltrami (LB) operator, and the inverse capacitance matrix is equivalent the Green function 
of this differential operator in the appropriate continuum limit. Let us denote this Green 
function by C _1 (x, x') and define U (x, x') = limc^o C _1 (x, x'). In a flat space C _1 (x, x') = 



K (2eyCi/Co |x — x'|) where Kq(x) is the zero order Bessel function of imaginary argument, 
and C/(x,x') = (27rC 1 )- 1 log|x-x'|.lli 

If we take the continuum limit where a becomes an infinitesimal then we can replace the 
sums by integrals over area and write the action as: 



1 rP r 

S[n(x),<^(x)] = -- jf y dxg 



1/2/ 



, -in(x) (Wx) + 2ttc; (0) (x)) - 

n Jo J 1 v ' 

+ § (^(x) - A-(x) + 2tt4 ) (x)) (^(x) - AM) + 2tt4 ) (x)) g* 



2e 2 j chd ^ 1/2 (x) n(x) n(x') C'^x,^) 



(16) 



where n(x a ) is the charge density. 

We now wish to perform a series of operations on the action so as to bring the vortex 
degrees of freedom to the fore. These steps are standard ones for deriving the interactions 
of vortices in 2D systems, made slightly more tedious by the necessity to keep track of the 
local metric. While they appear messy and involved, one should keep in mind that these 
steps are identical to what one would do if we had a regular, uniform array, but wished to 
evaluate all quantities in curvilinear (i.e. polar) coordinates. Furthermore, if we choose our 
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lattice to be smooth and undistorted, so that gij = 8ij, we should regain the standard vortex 
interaction. 

We first decouple the Josephson term via a Hubbard- Stratonavich transformation, intro- 
ducing a new field p(x). The second term then becomes: 

1 



IE, 



PiPj9 jk + {d k f ~ A k + 27rv^) Pj g jk 



(17) 



The action is now linear in ip; integrating it out we obtain the constraint: 



n(x) = -g- 1 ' 2 d k g x / 2 g* Pj 



(18) 



The right hand side is just the divergence of p expressed in invariant form. In order to 
satisfy this constraint we first express it in a "space-time" notation, using greek indices 
ji G {0, 1,2}, setting do = d t and defining the 3-momentum = (n(x),pi(x),p 2 (x)). Then 
eq.(|18D can be rewritten 



sT D^Pu = o 



(19) 



where g^ v is the the 3x3 metric, 



\ 9ij J 

and is the covariant derivative. Note that trivially detg^ 
satisfied if p M is the covariant curl of a vector field: 

e^g^DxK^ = 2n Pu 



(20) 



det gij. Eq.(|T9|) can be 



(21) 



where is an auxiliary or gauge field and e KAAt is the totally antisymmetric Levi-Civita 
tensor. Since covariant derivatives of the metric are zero, eq.(plj) is satisfied automatically. 
We rewrite our action in terms of the gauge fields K^, expressing n(x) and p(x) using eq.(^Il) 
As is common with such fields, we may choose a gauge condition for that will simplify 
the algebra: we choose 



S 



g i W i K j = 



(22) 



in which case we the action becomes: 



S[K^} = --jdrjd^ g l ' 2 {x) [ie^ x {v x D^K v - A x D^K U 



1 



OTX Cj J 



'J 



- J d^g x l\^)e %] e u A^(x) D k K t (x!) CT^x') 



2tt 2 C 



where = (0, A%, A 2 ). Integrating by parts on the first term we obtain 



(23) 



=-i I 'It /r/xV 2 (x) 



i<T (J v - $„) ^ + 7t4- (5^o djK + doK, d Kj] 
ihjj 



+ J cbdgWWeijektDiKjW D k K t {x!) C~\^) 



(24) 



where $ is a magnetic field density, $ M = (S z /$o,0,0), and J„ is a vortex current density 
related to our field v by a covariant curl:E3 



(25) 



This is the continuum version of the distorted Josephson junction lattice. The field 
mediates the interation between vortex charges and currents. The action is quadratic in 
these fields, and thus the path integrals over can be performed. The resulting action will 
give the vortex Hamiltonian in the curved space. This is done in section [TTJ below. 

III. VORTEX INTERACTIONS IN A CURVED SPACE 



The continuum action of eq. (|24| ) is expressed in terms of the gauge fields and the vortices. 
We would like to integrate out the so as to have an action solely in terms of the vortex 
degrees of freedom. First, we write the vortex current density as a sum over discrete point 
vortices located at positions X^ n ^: 



^ = Ep( x<n) W) tf(x-x<»>(t)) (i,xi n \t),xf\t) 



(26) 
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where p(X( n )) is the charge of the vortex located at X^ n \ 

We can integrate out the K field as follows: integration by parts on the quadratic 
term yields K§g~ x l 2 dig x l 2 g %3 dj Kq\ this combination of metrics and derivatives is just the LB 
operator, the generalization of the Laplacian to this curved space. Integrating out Kq will 
yield a potential interaction between vortices, Ufa, x'). 

We next examine the dependence. The d t Ki terms correspond to a retarded interac- 
tion between the vortices; we can neglect this non-locality in time so long as the characteristic 
velocity of vortices is much less than uiji, where ujj = y/8EjEc/h, and i is a characteristic 
island size.0 Doing so amounts to making the charge and vortex degrees of freedom exactly 
dual. We further assume that Cq is small and we can approximate the inverse capacitance 
kernel with Ufa, x'). As discussed above, Ufa, x') is the inverse of the LB operator; when 
we perform the path integral it must cancel out the powers of V acting on Ki. This inte- 
gration will generate two types of terms quadratic in the vortex velocities. The first will be 
self-interaction terms that generate a mass for the vortex; the second terms are a velocity 
dependent interaction between the vortices.0 Our final action then depends only on the 
vortex degrees of freedom 

Putting this all together we get: 

J ^ n ^ m,n 

+- J2 p(X (m) ) p(X(")) U(p rn \ XW) (27) 

^ m,n 

where m c g = ix 2 h 2 Ci/2e 2 and we have suppressed the explicit dependence on time. The first 
term is the kinetic energy of the vortex in the twisted space; the second term is the veloctity 
dependent vortex-vortex interaction, and the third is vortex-vortex potential. 

In order to demonstrate the effect of the defects, we calculate the interaction kernel 
Ufa, x') for two simple cases, that of the disclination and the dislocation. In order to do 
so, we are simply solving for how point charges interact in the 2D curved space defined by 
the distorted lattice. In flat space the field about a point charge falls off logarithmically. If 
we set z\ = x + i y, the 2D potential of a point charge located at Z in a flat plane can be 
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written as: 

V(x,y) = V(z) = ^-Req\n(z-Z). (28) 

We will need this formulation below. The interaction energy for a set of charges in flat space 
is then: 

£flat = ^£p(X«) p(XW) Iog|X« -XW| (29) 

hJ 

In this result we have neglected a constant self-energy term for the vortices that diverges as 
the log of the system size. This term either forces the net vorticity of the system to be zero, 
or is compensated by the externally applied field. 

Although this is a continuum model for a problem on a lattice, it is known that this 
Green function is well approximated by that of the continuum problem!!! This gives us 
confidence in similar appoximations we make when we examine lattices with defects. 



A. Disclinations: 

A negative disclination (fig(|l|b)) in 2D can be viewed as the projection of a cone into the 
plane. In order to form the cone, a wedge of sites is removed from the lattice. Next, the 
exposed edges are brought together, puckering the surface into a cone. Finally, the lattice is 
projected from the 3D cone back down to a 2D plane. This final distortion or stretching of 
the lattice is irrelevant from the point of view of the excitations, since it does not change the 
number or connectivity of the junctions. We can thus approximate the interaction energy 
of a charge or vortex excitation with a negative disclination by calculating their energy on a 
cone, or equivalently, on a plane with a wedge removed (fig. 0). We therefore need only solve 
the Poisson equation for a single point charge on a cone in order to obtain the potential 
f/(x,x'). 

This problem on a cone is equivalent to a 2D electrostatics problem of a single point 
charge on the planar surface with a missing wedge, as depicted in fig. ^. The boundary 
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condition across the cut is that the electric field is continuous. We are free to choose our cut 
in any direction; we choose to place our charge on the positive x-axis and center the cut out 
wedge along the negative x-axis. The boundary condition now simplifies to the requirement 
that the normal component of the field is zero on the edges of the wedge, and the parallel 
component is continuous. We can achieve this by a simple conformal mapping.0 We can 
open a wedge in the plane via the mapping z — > z' = z x l p For a wedge with an opening 
angle a (fig. ||D we require p = 27r/ (27r — a). Then 

V{z) = Re \n(z p - Z p ) (30) 

From this result we can obtain the vortex-vortex interaction energy. In addition, the self- 
energy term generates an interaction energy between the charge and the disclination. This 
can be found by calculating the electric field near the charge, subtracting the part due to 
the self-interaction, and integrating the force on the charge to get an effective potential. 
Putting the two together we find that the interaction energy of a set of vortices (or charges) 
on a lattice with a disclination at the origin is: 

£disc = ^ IE P( xW ) P( x0) < 

log [\X (i) \ 2p + \X (j) \ 2p - 2|X W | P |X (J) | P cosp( 7i - 7i )" 
£p(X«) 2 (p-l) lnX«| (31) 

with 7 = arctan(X 2 /X 1 ). Eq. [H] reduces to the correct result for p = 1 (no distortion) and 
p = 2 (a half plane with an image charge at —z ). Translation invariance is explicitly broken 
by the defect. 

We can see a intuitively why such the electrostatic charge interacts with the defects by 
considering the density of field lines. Imagine a series of concentric circles drawn about the 
point charge. If a given circle does not intersect the wedge, then the circumference of the 
circle of radius R will be 2irR. Once the circle intersects the tip of the wedge (or the apex 
of the cone, if we think in 3D), the circumference of the circle will be smaller than 27rR, and 
the density of field lines radiating outward will be slightly higher. Thus the electrostatic 
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energy of the charge is lowered by moving in farther away from the tip of the cone. The 
reverse applies to a positive disclination, which allows the field to less dense, so that the 
electrostatic charge is attracted to the defect. This effect will not affect the vortex- vortex 
(or charge-charge) interaction if the particles are much closer to each other than they are to 
the disclination, ( |X« - X & | < |X W | , | X & \ 

The vortex- vortex interactions are altered by the presence of the defect, introducing 
terms where the topological "charge" is coupled to two such excitations. However, we do 
not find that three body interactions are introduced among the vortices themselves. 



B. Dislocations: 

A similar calculation can be performed for a dislocation of Burger's vector b which we 
choose to be located at the origin and corresponding to an additional row of junctions along 
the negative x-axis, b = b e.2 as in fig.(|l}a). In this case it is easiest to start with the 
solution for an infinite half-plane, setting p = 2 in eq(|30D. We then deform the half plane 
into an semi-infinite slot of width b extending along the positive x-axis using the general 
Schwarz transformation.0 The resulting exact mapping is not invertible in a closed form; 
we approximate it with: 

b 

z -> z + —\nz (32) 
2vr V ; 

where we take | argz| < 7r. This result is asymptotically correct for |b| <C |r|, but lacks the 
sharp corners of the exact transformation. In this limit we find 



£ di si^Ep(x (i) ) p(x (j) ) 1^ 
-Ep(x (i) ) 2 ^ 



Xf > - X? - A( 7 W - ^))\ 2 + ( X« - + - In 



2tt v ' ' ') \ 1 1 2vr Xb 



to lowest order in b. In this expression X = X\ ei + X2 62, and the dislocation is at the 
origin. The interaction energy between a charge and a dislocation at the the point x is in 
general 
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to lowest order in b. We see that vortices are attracted or repelled from a dislocation as with 
a disclination, but the effect is weaker and depends on the orientation of the dislocation as 
well as its distance. The effect is smallest when b is perpendicular to the field lines. 

A metric with such a dislocation is said to possess torsion. Such metrics appear in general 
relativity with spinning masses.0 



IV. MOTION OF VORTICES IN A CURVED SPACE 

We now examine the actual motion of a single vortex in the above two cases. Our problem 
simplifies considerably because we do not have to worry about vortex-vortex interactions. 
In order to proceed we need to specify the metric. For the disclination, we choose to work 
in polar coordinates. The metric can be written 

d o ^ 

(1 



9ij 



\ 



QL\2 2 
2n> ' 



(35) 



For a = (or p = 1) we recover the flat space metric. Note that the circumference of a 
circle of radius R in this space is (2n — a)R. Our Lagrangian in polar coordinates (r, 7) is 
then 



L = \ m cff gij Xi Xj - ^ In A 

m eS ( . 2 . /1 a \2 2-2 



r 2 + (1 - 2^7' 



— (p— 1) lnr 



(36) 



Solutions of this system can be reduced to quadrature: 

(7 - 7o) 



lnr 



(37) 



where C\ and C2 are constants defined by the initial energy and angular momentum. In fig.(f§) 
we show some sample trajectories for vortices propagating ballistically in a junction array 
with a disclination. Note that for positive disclinations we see a repulsion, and for negative 
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ones we see an attraction. Also shown are trajectories for fictitious "neutral" vortices that 
propagate along straightlines. These trajectories show the effect of the vortex-disclination 
potential. 

For the dislocation with Burger's vector b = b e\, located at the origin, we choose to 
work in Cartesian coordinates. We approximate the metric to lowest order in b: 

( , by 



9ij 



by 
V 2^7 



2ttt 2 
, 2bx 
2nr 2 



) 



(38) 



In polar coordinates the Lagrangian for the vortex is 
1 



r + r 7 + 2by (f sin 7r7 cos 7) 



bVo 
2tt 



cos 7 



(39) 



For b = 1 we recover the flat space metric. In fig.® we show some sample trajectories 
for vortices propagating ballistically in a junction array with a dislocation. The trajectory 
depends upon the relative orientation of the dislocation and the vortex. Also shown are 
trajectories for "neutral" vortices that propagate along straightlines. These trajectories 
show the effect of the vortex-dislocation potential. 



V. CONCLUSION 

In this paper we have shown that the geometry or connectivity of a lattice can have 
a strong affect on the dynamics of vortices within the lattice. These effects can be seen 
in the ballistic propagation of vortices. Such experiments have already been performed on 
regular lattices. We propose that similar experiments could be done with lattices possessing 
topological defects. 

In addition, defects might also change the thermodynamic properties of the lattice. The 
presence of a postitive disclination acts as an attractive potential for both positive and nega- 
tive vortices, and it alters the interaction between vortices whose separation is comparable to 
their distance to the defect. These effects may both alter the Kosterlitz-Thouless transition 
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within the lattice. In addition, it will alter the interactions in the quantum Hall analogs 
proposed in such arrays.llM! 

Since the charge and vortex degrees of freedom are approximately dual,00 we expect 
similar results for charge excitations within the lattice. In the limit that we neglect the self- 
capacitance and the inertial terms in K t the two exactly dual. The interaction potentials 
for charge excitations will be identical to those calculated above, since they are both based 
on finding the inverse of the LB operator. Thus we expect that charge excitations on the 
lattice will also interact with distortions of the lattice. However, the Coulomb interaction 
is not strictly two dimensional in real lattices. We expect therefore that interactions with 
lattice defects will effectively be screened out at distance that depends upon the ratio of the 
nearest neighbor capacitances to the next-nearest neighbor capacitances. 

The motion of charges on twisted lattices also has potential application in tight-binding 
models of lattices with dislocations, or of carbon nanotubes. In these cases the lattice gives 
rise to metrics for the kinetic energy that have non-trivial torsion and curvature. These may 
act as sources of additional scattering or resistance. 

Similar interactions between defects in different "order parameters" has been proposed 
in other systems. This idea is at the heart of the superhexatic stated! in which the defects 
that mediate the melting transition interact with superfluid vortex degrees of freedom. A 
mathematically similar situation arises in the case of tilted smectic liquid crystals. @ This 
system is composed of a 2D surface of rod-like molecules that orient themselves at a small tilt 
with respect to the surface normal. The molecules themselves form a liquid crystal, which 
has disclinations as spatial defects. The U(l) symmetry associated with the orientation of 
the tilt interacts with the spatial defects in the smectic. The interaction energy between the 
two (sometimes called "white" and "green" vortices) is analogous to the vortex-disclination 
energy we have discussed above. It is also amusing to note that these Josephson junction 
lattices also provide analogs of how charges interact in a 2D gravitational system. 

It is quite common to model continuous differential equations by putting them on a mesh. 
It is usually assumed that the nature of a mesh has no effect on the dynamics of the system, 
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if the mesh is sufficiently fine. We have shown that topological properties of the mesh can 
introduce effects in the dynamics of the system being modelled. Care must be taken when 
using such non-integrable metrics, (i.e. those with torsion or curvature). 

In the world of superconducting lattices, distance is measured in junctions. This allows 
the experimenter to create lattices where the excitations move in a curved space, producing 
several novel effects. 
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FIGURES 

FIG. 1. (a) A square lattice with a dislocation. It can be viewed as a missing horizontal row of 
islands along the positive x direction, or an inserted row of island along the negative x-axis. (b) A 
square lattice with a negative disclination. It is formed by removing a 90° wedge from the original 
square array, and then distorting the lattice so that the open edges meet. Although the system has 
a center of three-fold symmetry, note that all sites still have four neighbors. To excitations within 
the lattice it will appear to be a curved space. A positive disclination (not shown) is formed by 
the insertion of a 90° wedge of material into the lattice. It is also possible to form disclinations in 
triagular lattices by removing or inserting a 60° wedge of material. 

FIG. 2. A schematic representation of a junction lattice where each island is treated as a point. 
With each lattice point we associate two vectors pointing to the nearest neighbors to the left and 
above the island. The vectors u^) and U( 2 ) may vary as we move through the lattice. 

FIG. 3. To excitations moving within the a junction array, a lattice with a disclination is 
equivalent to the surface of a cone. If we slice the cone from the tip out to infinity, we can flatten 
it to the plane. Here we sketch the field lines of a point charge P in such a case, where the wedge 
of removed material forming the disclination has an opening angle a, and is centered at the origin, 
O (the tip of the cone). Derivatives of the field are continuous across the cut so by symmetry, 
the field will approach it tangentially. Because the field lines are "excluded" from the wedge (the 
shaded area), the field lines are forced closer together, increasing the energy. This results in a 
repulsive interaction between the point charge and the defect. 

FIG. 4. Sample trajectories for a vortex moving in a disclinated lattice, as in fig.(l.b). The 
disclination is at the origin and has a "charge" or missing angle of 7r/4. Show here are trajectories for 
vortices with an initial velocity vq and an impact parameter of r^. The dashed line is a trajectory 
for a "neutral" vortices showing just the effect of geometry. The other lines are trajectories of 
vortices with interaction energies T = Ej/m c gVQ of 1/4, 1/2, and 1, from top to bottom. The 
vortices start at the position (ro, — vq). The distances are measured in tq. 
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FIG. 5. Sample trajectories for a vortex moving in a lattice with a dislocation, as in fig.(l.b). 
The dislocation is a missing row of junctions along the positive x-axis, and is centered at the origin. 
Show here are trajectories for vortices with an initial velocity vo and an impact parameter of ro. 
Plotted are trajectores for "neutral" vortices (dashed line), and vortices with interaction energies 
T = Ej/m e fiVQ of 0.40 (middle line), and 1.0 (outer line). The distances are measured in ro, and 
b/r = 0.10. 
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Figure la 
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Figure lb 
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Fig. 3 
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